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SUMMARY: 

Modern liquid rocket engines often use liquid oxygen (LOX) and liquid hydrogen (LH 2 ) as 
propellants to achieve high performance, with the engine operational conditions in the supercritcal regimes 
of the propellants. Once the propellant exceeds its critical state, it essentially becomes a puff of dense fluid. 
The entire field becomes a continuous medium, and no distinct interfacial boundary between the liquid and 
gas exists. Although several studies have been undertaken to investigate the supercritical droplet behavior 
at quiescent conditions, very little effort has been made to address the fundamental mechanisms associated 
with LOX droplet vaporization in a supercritical, forced convective environment. The purpose of this work 
is to establish a theoretical framework within which supercritical droplet dynamics and vaporization can be 
studied systematically by means of an efficient and robust numerical algorithm. 

TECHNICAL DISCUSSION: 

Several analyses based on the model of a spherical droplet in a quiescent environment were recently 
developed to address the characteristics of supercritical droplet vaporization([l][2][3][4]). The model is 
adequate only for small droplets which become equilibrium with the ambient flow. Contrarily, convective 
effect is clearly observed for large droplets as a result of high momentum inertia of the liquid phase. Figure 
1 depicts the situation examined here, showing an isolated vaporizing LOX droplet in a convective hydrogen 
environment. The initial temperature of the droplet is subcritical, but the ambient temperature and pressure 
are in the supercritical regime of oxygen. As a result of heat transfer from the surrounding gases, the droplet 
surface temperature increases rapidly and reaches its critical mixing point. Under this condition, the sharp 
distinction between gas and liquid phases disappears. The droplet regression can be best characterized by 
tracing the location of the critical temperature. 


Governing Equation 


The analytical model is based on the complete time-dependent conservation equations of mass, 
momentum, energy and species concentration for a multicomponent system with variable properties. If 
body forces, viscous dissipation, radiation, and chemical reaction are ignored, these equation can be written 


in the following vector form. 


Mass: 


Momentum: 


Energy: 


Mh + Jh iA ‘= 

( 1 ) 

ijjjpu^V + JJpuiUjdAj = JJnjdAj. 

( 2 ) 

JJJp edV + JJp eu i dA i = JJ T ij u idAj - JJ(q T )jdAj. 

( 3 ) 
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Species Concentration: 


where 


Jt IJJ pYldV + JJ pY,u > dA i = ~ . (4) 


Tij = -pSij + 2 neij - -/i(V • u)6 i} . 

Standard notations in fluid mechanics and thermodynamics are used in (l)-(4). The specific total energy e 
is defined as 


1 t T 

1 = 1 JT re 


V ti* tii 

c '^-? + -r 


where the index / represents the number of species considered, Y/ the mass fraction of species /, and T re j the 
reference temperature for energy. Fick’s and Fourier’s laws are used to approximate the species and thermal 
diffusion in (3) find (4), respectively. 


Property Evaluation 

Thermophysical properties, including both thermodynamic and transport properties, play a decisive 
role in determining droplet vaporization behavior. It is generally accepted that thermal conductivities and 
heat capacities of pure fluids and fluid mixtures can be divided into three contributions, and correlated in 
terms of density and temperature. For example, the thermal conductivity of fluid oxygen can be expressed 
as [5] 

A (p,T) = A 0 (T) 4- AA txe (p,T) + A\ crit (p,T). (5) 

The first term on the right-hand side represents the value in the dilute-gas or zero-density limit, which is 
independent of density and can be accurately predicted by kinetic-theory equations and an m-6-8 model 
potential. The second is the excess thermal conductivity which, with the exclusion of unusual variations 
near the critical point, characterizes the deviation from A 0 for a dense fluid. The sum of A 0 and AA CXC is 
sometimes called the “background” thermal conductivity. The third term refers to the critical enhancement 
which accounts for the anomalous increase above the background thermal conductivity as the critical point 
is approached. There fire two separate regions where the critical enhancement of thermal properties is 
important. The first region is located in the close vicinity of the critical point, and is the region where 
a scaled equation of state must be used, as defined by Sengers et al[6]. For oxygen with T c = 154.581 K 
and p c = 0.436 g/cm 3 , this region is bounded approximately by 150 < T < 160 K and 0.32 < p < 0.544 
g/cm 3 . The second region, designated as the extended critical region, covers the p-T domain in which critical 
enhancement is significant, but the fluid behavior can be described by an analytical equation of state such 
as a modified Benedict- Webb- Rubin (BWR) type. This region can extend to quite high temperatures, up to 
approximately 2T C . The critical enhancement along various isotherms can be centered on a density, p C ente n 
which deviates slightly from the critical density and decreases with increasing temperature. More detail 
description can be found in [4]. 


Numerical Algorithm 

The model considered here is characterized by strong couplings among heat transfer, transient 
diffusion, and phase transition, etc.. Because various time and length scales are involved in the entire process, 
the numerical scheme for solving this serious stiffness problem becomes a challenge. Moreover, the singular 
behavior of momentum balance in the limit of diffusion-dominant process adds further complications. In view 
of these difficulties, a numerical algorithm capable of treating time-accurate low Mach-number compressible 
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flow has been developed. The scheme is constructed in two steps. First, a rescaled pressure term is used in 
the momentum equation to circumvent the pressure singularity at very low Mach number. 

p(x,y,t) = p 0 +p(x,y,t) 

The second step is based on the establishment of a dual time-stepping integration procedure. With the 
addition of a set of well-conditioned artificial time derivatives to the conservation laws, a new system of 
governing equations is obtained whose converged solution in pseudo-time corresponds to a time-accurate 
solution in the physical- time domain. 

r dZ dQ djE-K) d(F — K) 

dr dt 8x + 8y 

where r presents pseudo time, and Z = (p', u, v ) h t Yi). 

RESULTS: 

As a first study of supercritical droplet behavior, A spherical droplet with an introductory 
temperature of 300 K is first placed in a supercritical, forced convective environment (Too = 1000 K). The 
droplet Reynolds number based on the relative velocity is 60. The droplet is 100 pm in diameter and 
contains 100% oxygen initially. As a result of large temperature gradient between gas and liquid phases, 
surface temperature raises and reaches critical condition shortly. Due to the dimishment of surface tension, 
droplet deformation and shattering can be clearly observed. Figure 2 shows the density contours, streamlines, 
and velocity vectors at t = 0.02ms in a supercritical environment with Hi as the ambient gas. Figure 3 
shows the oxygen concentration contours. Results indicate that the ambient gas mass diffusivity plays an 
important role in determining the droplet behavior. A high mass diffusivity causes rapid penetration of the 
ambient gas into the droplet, and consequently modifies the droplet dynamics (e.g., deformation, shattering, 
etc.) caused by the density difference between the droplet and the gas. 
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Figure 1: Schematic of a vaporizing LOX droplet in hydrogen stream. 



Figure 2: Density contours, streamlines and velocity vectors of LOX droplet vaporizing in a hydrogen 
environment. Re = 60, t = 0.02ms. 


31 


Mass Diffusion Exerts Strong Influence on Droplet Behavior in 
a Supercritical, Forced-Convective Environment 
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Figure 3: Time history of mass fraction of a vaporizing droplet at supercritical condition. 



